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The shell model in the complex /c-plane (the so-called Gamow Shell Model) has recently been 
formulated and applied to structure of weakly bound, neutron-rich nuclei. The completeness rela- 
tions of Newton and Berggren, which apply to the neutron strictly valid for finite-range 



potentials. However, for long-range potentials, such as the Coulomb potential for protons, for 



which the arguments based on the Mittag-Leffler theory do not hold, the completeness still needs 

h ; 

to be demonstrated. This has been done in this paper, both analytically and numerically. The 
generalized Berggren relations are then used in the first Gamow Shell Model study of nuclei having 
both valence neutrons and protons, namely the lithium chain. The single-particle basis used is that 
of the Hartree-Fock-inspired potential generated by a finite-range residual interaction. The effect 
of isospin mixing in excited unbound states is discussed. 
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I. INTRODUCTION 



One of the main frontiers of the nuclear many-body problem is the structure of exotic, 
short-lived nuclei with extreme neutron-to-proton ratios. Apart from intrinsic nuclear struc- 
ture interest, properties of these nuclei are crucial for our understanding of astrophysical 
processes responsible for cooking of elements in stars. From a theoretical point of view, 
the major challenge is to achieve a consistent picture of structure and reaction aspects of 
weakly bound and unbound nuclei, which requires an accurate description of the particle 
continuum Q|. Here, the tool of choice is the continuum shell model (see Ref. Q] for a 
recent review) and, most recently, the Gamow Shell Model (GSM) 0, U, 0, Q (see also 
Refs. 

ELBE]). 

GSM is the multi-configurational shell model with a single-particle (s.p.) 



basis given by the Berggren ensemble |l( 
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which consists of Gamow (or resonant) 



states and the complex non-resonant continuum. The resonant states are the generalized 
eigenstates of the time-independent Schrodinger equation which are regular at the origin 
and satisfy purely outgoing boundary conditions. The s.p. Berggren basis is generated by 
a finite-depth potential, and the many-body states are obtained in shell-model calculations 
as the linear combination of Slater determinants spanned by resonant and non-resonant s.p. 
basis states. Hence, both continuum effects and correlations between nucleons are taken into 
account simultaneously. The interested reader can find all details of the formalism in Ref. 
^J, in which the GSM was applied to many-neutron configurations in neutron- rich helium 
and oxygen isotopes. 

When extending the GSM formalism to the general neutron-proton case, with both pro- 
tons and neutrons occupying valence s.p. states, one is confronted with a theoretical problem: 
the Berggren completeness relations, which are the pillars of the Gamow Shell Model, have 
been strictly proved (and checked numerically ^, ^|) only for quickly vanishing (finite- 
range) local potentials, while the repulsive Coulomb potential for protons has infinite range. 
The theoretical problem lies, in fact, not in the Berggren (complex-energy) completeness 
relation itself, but in the Newton (real-energy) completeness relation jl4 Jig. This latter 
involves both bound and scattering states, upon which the Berggren completeness relations 
can be demonstrated using the method of analytic continuation. 

Our paper is organized as follows. Section |H] contains a derivation of the Newton com- 
pleteness relation that is valid for a rather wide class of potentials, including the Coulomb 
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potential. Based on this result, the Berggren completeness relation for protons is derived in 
the same way as previously done for neutrons [4]. The numerical tests of the completeness 
of the proton Berggren ensemble is given in Sec. 11111 Section HV1 introduces the Hartree-Fock 
(HF) inspired procedure used to optimize the s.p. basis, the so-called Gamow-HF method. 
The first GSM calculation involving active neutrons and protons is presented in Sec.|VJ with 
the the lp-shell study of the lithium chain, ranging from 5 Li to 11 Li. The residual interaction 
used is a surface-peaked finite range force. A novel aspect, absent in our previous GSM stud- 
ies, is the appearance of T=0 couplings which seem to exhibit significant particle-number 
(or density) dependence. Finally, Sec. I VII contains the main conclusions of our work. 

II. COMPLETENESS RELATIONS IN THE GSM: ANALYTICAL CONSIDERA- 
TIONS 

As the one-body completeness relation for resonant and scattering states is prerequisite 
for our theory, we shall demonstrate it rigorously. We shall first consider the case of a local 
potential and then generalize it to nonlocal potentials. 

A. Local potential 

In order to demonstrate the orthonormality and completeness relations for s.p. proton 
states, we consider a spherical proton potential that is finite at r=0, and it has a pure 
Coulomb behavior for r — > +00. The one-body radial wave functions u(r) are solutions of 
the Schrodinger equation: 



where potential v is given in units of fm -2 , and I is the angular momentum of the particle. 
Let us consider the bounded region enclosed in a large sphere of radius R. (This can be 
achieved by introducing an infinite well of radius R surrounding the nucleus.) Of course, in 
the final result, R will be allowed to go to infinity. For each value of R, one has the following 




(1) 



(2) 



r 
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completeness relation on [0 : R] 



[16]: 



+00 



J> n (r> n (r') + Yl ^ d \k m ,r)u {d \k m ,r') = S(r - r'), (3) 

neb m=0 

where b denotes the set of bound states having radial wave functions u n {r) with k\ < v(R), 
and u^(k m ,r) is a wave function of a normalized discretized continuum state, given by the 
boundary conditions u^(k m , 0) = and u^ d \k m , R) = 0. 

For the purpose of this discussion, it is convenient to introduce the set of wave functions 
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. u {d \k h r) 

u(k u r) = -7===, (4) 
V ki+i — ki 



which obey the following normalization condition: 



(u{ki)\u(kj)) = -j—^—jT- 

Ki+l Ki 



Since, in addition, 



(u(k)\u n ) = (6) 
(u n \u n >) = 8 nn >, (7) 

the box completeness relation can be written as: 

+00 

^u n (r)w n (r') + ^ u ( k ™, r)u(k m , r')(k m+1 - k m ) =5(r-r'). (8) 

neb m=0 

When R — > +00, the infinite series in Eq. (JHJ) becomes an integral, thus giving the expected 
completeness relation. Unfortunately, this cannot be done right away, as the series and the 
integral converge only in a weak sense. 

To prove the convergence rigorously, let us consider the completeness relation of the free 
box expressed in the form of Eq. 

+00 

^ fii?i(v)j((v')(Km + i - K m ) = S(r - r'), (9) 

where ji(n m r') are the Riccati-Bessel functions, ji(n m R)=0 (rn— 0,1,2...), and B m is the 
normalization constant: 

-bI^I+^rYj^^r) = — l - — . (10) 



Subtracting Q from (JHJ), one obtains: 



y]u n (r)u n (r') 

+00 

m , r)tt(/c m , r ) (fc m +i fc m 

) - B mjl( K mr)jl( 



m=0 

0. 



(11) 



We shall now demonstrate that the above series converges in the sense of functions, so the 
limiting transition from a series to an integral when R —>■ +00 can be easily carried out. 

To this end, let us consider the behavior of the m-th term in the series when m (and k m ) 
— > +00. For very large values of k m , one can use the semiclassical expansion in powers of 
kZ 1 : 



^(kmj r ) — C m ji(k m r) — C rn — — ji{k m r) + O ( ~r^~ 



ji{k m r) 



71 



sin [ k m r - I— 



a i I , n \ ^ I 1 
cos k m r — /— + U 



*2km.v 



kl 



2 / ' 



(12) 
(13) 



where C m is a normalization constant, a\ is a constant depending on I only, and 

rr 

V(r) = / f(r') dr'. (14) 

For the Coulomb potential, V(r) oc lnr; hence the expression (II 2|) properly accounts for the 
logarithmic term in the phase shift. It immediately follows from Eqs. (|12p and (|13p that 



/-„. ; (m + _j)> + „ + *V(*) +0 fJ^ 



i2 2^77777 

(m + |) 7r a ; 



+ 



i? 2Rm7T \m £ 

The constant Cm can be determined from the normalization condition: 



(15) 
(16) 



Cl 



r r 



j\ (k m r) 



V(r] 

h 

A, n, 



ji(k m r)j'i(k m r) 



C 2 

' h ()i lf 



b m+l 



(17) 



Since the integral involving V behaves like 1/k^, C m becomes: 
R 



j l (k m R) - ji+i{.k m R)ji-i{.k m R) 



h , — h ■ 



(18) 



cf. Eq. (Unj). 
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Using Eqs. (|T5j) and (JTHJ), one obtains: 



Cm = \ll + 0(l-\, (19) 



7r \m 2 

*" = v§ +0 Gy- < 20 » 

Let us consider the behavior of B m for — > +00 but K m — > k with k>0. The expansion of 
Eq. ()13|) is still valid in this for r=R, it is a familiar expansion in 1/R of the Bessel 

function. It follows from Eq. (fT3*|) that: 



K1K ofi| (21) 



i2 V^ 2 

with m chosen SO is the closest to k. Then, Eqs. (jl(Jj) and (|13j) give: 

s -=y! +o (i)- < 22) 

as expected. 

Note that the leading term in Eqs. (|19|) and (|2(J|) is the familiar normalization of continuum 
wave functions The reminders, of the order of m 2 , guarantee the convergence of the 
series. By using Eqs. (fT2|) -([TT) j) . one can show that the series (|TT|) converges for all r > and 
r' > 0. 

As a consequence, in the limit of R — > +00, Eq. (jTTj) becomes: 



y2,u n {r)u n {r') + / 



2 ~ 

u(k,r)u(k,r') ji{kr)ji(kr') 

71 



dk = 0. (23) 



By taking advantage of the closure relation for the Riccati-Bessel functions, 

r + °° ~ ~ 71 

/ ji(kr) 3l (kr')dk=-5(r-r'), (24) 
Jo * 

one finally arrives at the sought completeness relation: 

POO 

^M n (r)«„(r') + / u(k,r)u(k,r') dk = 5{r — r'). (25) 
neb Jo 

By using the same arguments as in Ref . |4( , one obtains the generalized Berggren complete- 
ness relation, also valid for the proton case: 

J2 MfKH + f u(k,r H ky)dk = S(r-r'). (26) 

n£b,d J L+ 



For details, including the numerical treatment of scattering wave functions and correspond- 
ing matrix elements, we refer the reader to Ref. Let us only remark, in passing, that 
in the presence of the Coulomb potential the standard regularization procedure Q] has 
to be modified In our work, however, we apply the exterior complex scaling method 
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22j which works very well regardless of whether the Coulomb potential is used or not. 



B. Nonlocal potential 



In the presence of a nonlocal potential, such as the HF exchange potential generated by 
a finite-range two-body interaction, the Schrodinger equation (JTJ) becomes: 

'1(1 + 1) 



u"(r) 



+ Viir) — k 



u[r) + I v n i(r,r')u(r') dr' 
'o 



(27) 



where vi is the local part of the potential, and v n i its nonlocal kernel. We assume that 
v n i(r,r') — > when r — > +oo or r' — > +oo (nuclear potential has to be localized) and that 
v n i{f,0) = V r (the potential is regular at the origin). As the radial HF functions are 
regular at zero, the latter condition is automatically met for the HF exchange potential. 

If the integral containing the nonlocal potential v n i behaves like 1/k 2 when k — > +oo, 
then the asymptotic expression (fT^j) holds. Indeed, integration by parts yields: 

W)j'i{kr') 



'0 

1 

O 



ji(kr') 



^(r,0)Ji(0) 



k 2 



o 



2k 

12 



dr 



dr' 



r, r')Ji(kr') dr' 



+ 0| F 



(28) 



pr pt ft 1 

where Vj(r) = / vi(r') dr', Ji(t) = / ji(t") dt" dt', and to and t' are chosen so Ji(t) 

Jo J t Jt' 

is bounded on [0 : +oo[. Consequently, Eq. (|T2j) also holds for nonlocal potentials. The 

proof of completeness can be, therefore, performed in the same way as for local potentials, 

by simply replacing v by v\ in all expansions in k~ l . 
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III. COMPLETENESS OF THE ONE-BODY PROTON BERGGREN ENSEMBLE: 
NUMERICAL TESTS 



In this section, we shall discuss examples of t 
one-proton case (for the neutron case, see Ref. 



re Berggren completeness relation in the 
The s.p. basis is generated by the 



11C 

4] 



spherical Woods-Saxon (WS) - plus - Coulomb potential: 



V(r) = -V f(r) - W so l ■ s-^ + V c (r), (29) 

r dr 



f(r) 



1 + exp 



d 



(30) 



In all the examples of this section, the WS potential has the radius Rq=5.3 fm, diffuseness 
g?=0.65 fm, and the spin-orbit strength V so =5 MeV. The Coulomb potential V c is assumed 
to be generated by a uniformly charged sphere of radius i?o and charge Q=+20e. The depth 
of the central part is varied to simulate different situations. 

In this section, we shall expand the 2p 3 / 2 state, |wws) ; either weakly bound or resonant, 
in the basis |m ws b(/c)) generated by the WS potential of a different depth: 

\uws) = ^2c ki \u WS B(ki)) 

i 

+ / c(fc)|u WS B(fc)) dk, (31) 

JL+ 

cf. Eq. (J26|) . In the above equation, the first term in the expansion represents contributions 
from the resonant states while the second term is the non-resonant continuum contribution. 
Since the basis is properly normalized, the expansion amplitudes meet the condition: 

E c l+ / c 2 (k)dk = l. (32) 

In all cases considered, the 0p 3 / 2 and IP3/2 orbitals are well bound (by ~ 50 MeV and ~ 20 
MeV, respectively) and do not play any significant role in the expansion studied, although 
they are taken into account in the actual calculation. The 2^3/2 state is, however, either 
loosely bound or resonant, and the scattering states along the contour L + are essential to 
guarantee the completeness. To take the non- resonant continuum into account, we take the 
complex contour L + that corresponds to three straight segments in the complex k— plane, 
joining the points: k = 0.0 - i0.0, k x = 0.3 - iO.l, k 2 = 1.0 - i0.0, and k 3 = 2.0 - iO.O (all 



S 



in fm 1 ). The contour is discretized with n=40 points: 

n 

|2p 3/2 > ^ C «I M «) + ^2 c k>\ u ( k i))- ( 33 ) 

n£b,d i=l 

In the first example, we shall expand the 2p 3 / 2 s.p. resonance (£ =3.287 MeV, T =931 
keV) of a WS potential of the depth Vo=65 MeV in the basis generated by the WS potential 
of the depth V(P=70 MeV. (Here the 2p 3 / 2 s.p. resonance has an energy £ =1.905 MeV 
and width T =61.89 keV.) After diagonalization in the discretized basis (J33|) . one obtains 
E =3.289 MeV and T =934 keV for the 2p 3 / 2 s.p. resonance, i.e., the discretization error 
is ~3 keV. The density of the expansion amplitudes is shown in Fig. As both states are 
resonant, the squared amplitude of the 2p 3 / 2 basis state is close to one. Nevertheless, the 
contribution from the non-resonant continuum is essential. It is due to the fact that the 
resonant state in the basis is very narrow, whereas the expanded resonant state is fairly 
broad. It is interesting to notice that the contribution from scattering states with energies 
smaller than that of the resonant state is practically negligible; this is due to the confining 
effect of the Coulomb barrier. 

The second example, shown in Fig.El deals with the case of a 2p 3 / 2 state that is bound in 
both potentials. Here V =75 MeV and 1/ (B) =80 MeV, and the 2p 3/2 state lies at £=-0.0923 
MeV and £=-2.569 MeV, respectively. Here, the scattering component is almost negligible, 
which reflects the localized character of bound proton states. After the diagonalization, one 
obtains £=-2.568 MeV and T=1.73 keV for the 2p 3 / 2 state, which is indeed very close to 
the exact result. 

In the third example, the unbound 2p 3 / 2 state (E =1.905 MeV and T =61.89 keV; Vo=70 
MeV) is expanded in a WS basis containing the bound 2p 3/2 level (£=-0.0923 MeV; V K '=75 
MeV). As a consequence, the resonance's width has to be brought by the scattering states. 
Nevertheless, the component of the 2p 3 / 2 state of the basis is still close to one, whereas the 
continuum component plays a secondary role. Once again, one can see a Coulomb barrier 
effect: even if the expanded 2p 3 / 2 state is unbound, its wave function is very localized due 
to the large Coulomb barrier; hence it has a large overlap with the bound 2p 3 / 2 basis state. 
The diagonalization yields £ =1.9065 MeV and T =58.9 keV, i.e., the discretization error is 
again close to 3 keV. 

In the last example, shown in Fig. HJ the basis is that of the WS potential with a depth 
of 70 MeV, and the expanded state corresponds to a WS potential with Vo=75 MeV. This 
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■ ■ ■ imaginary part 
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1 
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FIG. 1: Distribution of the squared amplitudes c 2 (k) of the 2p 3 / 2 proton state of a Woods-Saxon 
potential with a depth Vo = 65 MeV, in the s.p. basis generated by a Woods-Saxon potential 
with a depth Vo = 70 MeV. The Coulomb potential is assumed to be that of a uniformly charged 
sphere. The amplitudes of both real (solid line) and imaginary (dotted line) components of the 
wave function are plotted as a function of $l[k]. The height of the arrow gives the squared amplitude 
of the 2p 3 /2 state contained in the basis. 

is the most interesting case since one expresses a bound (real) state in the basis which 
contains only complex wave functions. (The contribution from well bound and 
s.p. states is negligible.) As expected, the behavior of the amplitudes is very close to that 
of the previous example (see Fig. |3J). Moreover, one can notice that the scattering states 
become important in the expansion when their energies approach the resonant state energy. 
The diagonalization gives E =-0.0925 MeV and T =-3.8 keV for the 2p 3 / 2 state. 
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FIG. 2: Similar as in Fig. ^except for the bound 2p 3 / 2 s.p. state of the WS potential with Vo=75 
MeV expanded in the basis generated by another WS potential (Vq = 80 MeV). The height 
of the arrow gives the squared amplitude of the bound 2p 3/ / 2 state at the value of -9f[A;]. (The 
corresponding A;— value is purely imaginary.) 



Summarizing this section, our numerical tests demonstrate that the one-body Berggren 
completeness relation works very well in the proton case involving the Coulomb potential. 



For other numerical tests, see Refs. 



231 ] (study of the s.p. level density) and [24| (study of 



the Berggren expansion in the pole approximation). 
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FIG. 3: Similar as in Fig. [2] except for the 2p 3 / 2 resonance of the WS potential with Vq=70 MeV 
expanded in the basis generated by another WS potential (Vq = 75 MeV). 

IV. SPHERICAL GAMOW HARTREE-FOCK METHOD 

In our previous calculations of the He chain , we used the s.p. basis of the WS 
potential representing 5 He. This basis is appropriate at the beginning of the He chain, but 
when departing from the core nucleus, its quality deteriorates. For instance, the " 5 He" basis 
is not expected to be optimal when applied to the neutron-rich halo nucleus 8 He because of 
the very different asymptotic behavior of this weakly bound system. The obvious remedy 
is to use the s.p. basis that is optimal for a given nucleus, that is the HF basis. However, 
since the Berggren ensemble used in GSM is required to possess spherical symmetry, HF 
calculations must be constrained to spherical shapes. Moreover, since in some cases one 
is interested in unbound nuclei lying beyond the drip line (particle resonances), the HF 
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FIG. 4: Similar as in Fig. [2] except for the 2p 3 / 2 resonance of the WS potential with Vq=75 MeV 
expanded in the basis generated by another WS potential (Vq = 70 MeV). 

procedure has to be extended to unbound states. In the following, the HF-based procedure 
that meets the above criteria is referred to as the Gamow-Hartree-Fock (GHF) method. 

Since, strictly speaking, the spherical HF potential cannot be defined for open shell nuclei, 
one has to resort to approximations. In this work, we tried two different ways of averaging 
the HF potential. The first ansatz is the usual uniform-filling approximation in which HF 
occupations are averaged over individual spherical shells. In the second ansatz, the deformed 
HF potential corresponding to non-zero angular momentum projection is averaged over all 
the magnetic quantum numbers. Both methods reduce to the true HF potential in the case 
of closed-shell nuclei. 
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A. Average spherical HF potential 



In the uniform-filling approximation no individual HF orbitals are blocked. The matrix 
elements of the HF potential U u { between two spherical states a and j3 carrying quantum 
numbers (j, /) are: 

(a\U u{ \(3) = (a\h\[3) 

+ YTl S ^^(am Xm x \V\(3m Xm x ) (34) 

where h is the s.p. Hamiltonian (given by a WS+Coulomb potential), A is an occupied shell 
with angular quantum numbers l x ), N(X) is the number of nucleons occupying this shell, 
and V is the residual shell-model interaction. 

To define the M-potential, one occupies (blocks) the s.p. states in the valence shell that 
have the largest angular momentum projections on the third axis. The resulting Slater de- 
terminant corresponds to the angular momentum J—M. For closed-shell nuclei (M=0) and 
for nuclei with one particle (or hole) outside a closed subshell (M=j), this Slater determi- 
nant can be associated with the ground state of the s.p. Hamiltonian. However, in other 
cases it corresponds to an excited state with J>0. Spherical M-potential, Um, is defined by 
averaging the resulting HF potential over magnetic quantum number m: 

(a\U M \P) = (a\h\0) 

1 j 

+ — ^(am Xm x \V\(3m Xm x ) (35) 

l ' j m=j+l-N ltj x 

where Nij is the number of nucleons occupying the valence shell with quantum numbers 

The M-potential is expected to work better for nuclei with one particle (or hole) outside 
the closed subshell. However, one can expect this potential to be not as good as U u{ when 
the Slater determinant with J=M represents an excited state. 

B. Unbound HF states 

While the HF procedure is well defined for the bound states, it has to be modified for 
the unbound s.p. states (resonant or scattering), even in the case of closed-shell nuclei. 
First, the effective nuclear two-body interaction has to be quickly vanishing beyond certain 
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radius. Indeed, if it does not, the resulting HF potential diverges when r — > +00, thus 
providing incorrect s.p. asymptotics. Moreover, as resonant states are complex, the true 
self-consistent HF potential is complex. This is to be avoided, as the Berggren completeness 
relation assumes a real potential. However, since we are interested in the optimal basis- 
generating potential and not in the full-fledged complex-energy HF problem |25j, we simply 
take the real part of the (generally complex) HF potential. 



C. Treatment of the exchange part 

As the residual interaction used in our shell-model calculations is finite range, its exchange 
part gives rise to a non-local potential. The HF equations solved in the coordinate space 
can be written as integro-differential equations. The standard method to treat such a HF 
problem is by means of the equivalent local potential 26] : 

Veq(r) = vi{r) + -y ^ , (36) 

where we use the notation of Eq. ()27j) . The resulting HF equations are local but potentials 
become state dependent. The main difficulty is the appearance of singularities in V eq (r) 
due to the zeroes of the s.p. wave function. This problem is practically solved by replacing 
u(r) with a small number (e.g., 10~ 3 ) in the denominator of Eq. (J5fi)l when u(r) approaches 
zero, and by using splines to define the HF potential. The numerical accuracy is checked 
by calculating overlaps between different wave functions having the same (j, I) values. The 
overlaps are typically 10~ 5 , which is small enough to consider wave functions orthogonal. 

Let us note in passing that the demonstration of Sec. lIIBl is valid when the non-local part 
of the potential is localized, as it is the case for the nuclear interaction. However, if one wants 
to explicitly consider the Coulomb HF potential between valence protons, one has to resort 
to approximations in order to avoid its infinite-range non-local part. One possibility is to use 
the so-called Slater approximation, which has been shown to work fairly well [27 . Another 
method, the so-called generalized local approximation, has been proposed in Ref. 2j|, where 
the Coulomb exchange term has been parametrized in terms of a coordinate-dependent 
effective mass. In any case, the effect of the approximate treatment of the Coulomb exchange 
term on the GHF basis is very small as compared to other uncertainties related to the 
construction of the GHF Hamiltonian. 
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TABLE I: Binding energies of the He and Li isotopes (in MeV) calculated in the GSM using the 
GHF basis with (i) the uniform-filling approximation potential U u f or (ii) the LOvf-potential. See 
Sec. II V Al for definitions. 





6 He 


7 He 


9 He 


7 Li 


8 Li 


9 Li 


10 Li 


n Li 


U M 


-1.038 
-0.984 


-0.048 
-0.475 


-2.357 
-2.418 


-14.266 
-13.008 


-15.521 
-15.094 


-20.288 
-20.181 


-18.082 
-17.749 


-15.649 
-15.634 



D. Choice of the average potential 

In order to compare the quality of two GHF potentials, we inspect the binding energy 
(i.e., the expectation value of the GSM Hamiltonian) for several Li and He isotopes in the 
truncated GSM space. As discussed in Sec. |V| below, we took at most two particles in the 
GHF continuum. Due to this truncation, and as well as the discretization of the contour L + 
and the assumption of the momentum cut-off (the contour does not extend to infinity), the 
completeness relation of the resulting many-body shell-model basis is violated and, except 
for some special cases, the results obtained with different average potentials are different. 

According to the variational principle, better basis-generating potentials must yield lower 
binding energies. Table |l] shows the binding energies of 6 ' 7,9 He and 7_11 Li calculated in the 
GSM. Since 8 He and 10 He are closed-shell nuclei, both U n f and Um are identical in these 
cases. One can see that the uniform-filling approximation for the GHF potential works better 
for 7 Li, 8 Li, and 10 Li, whereas the s.p. basis of the M-potential is a better choice for 7 He. 
In all the remaining cases, the two potentials give results that are practically equivalent. In 
the particular case of 6 He and 6 Li (not displayed), there are at most two nucleons in the 
non-resonant continuum. Consequently, the shell-model spaces for 6 He and 6 Li are almost 
complete, and the results are almost independent on the choice of s.p. basis. Based on our 
tests, in our GHF calculations we use the uniform-filling approximation except for nuclei 
with one particle (hole) outside closed subshells where the M-potential provides a slightly 
more optimal s.p. basis. 
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V. GAMOW SHELL MODEL DESCRIPTION OF THE LITHIUM CHAIN 



A. Description of the calculation 

In our previous studies 0, Q|, we have used the s.p. basis generated by a WS potential 
which was adjusted to reproduce the s.p. energies in 5 He. This potential (" 5 He" parameter 
set 0]) is characterized by the radius R=2 fm, the diffuseness <i=0.65 fm, the strength 
of the central field Vo=47 MeV, and the spin-orbit strength V so =7.5 MeV. As a residual 
interaction, we took the Surface Delta Interaction (SDI). However, when it comes to practical 
applications, SDI has several disadvantages. Firstly, it has zero range, so an energy cutoff 
has to be introduced; hence the residual interaction depends explicitly on the model space. 
Moreover, as the SDI interaction cannot practically be used to generate the HF potential (it 
produces a nonrealistic mean field), one is bound to use the same WS basis for all nuclei of 
interest, which is far from optimal as the number of valence particles increases. So, we have 
decided to introduce p a finite-range residual interaction, the Surface Gaussian Interaction 
(SGI): 



V^\r 1 ,r 2 )=V (J,T)-ex V 



2 

ri - r 2 ■ 



•5(|r 1 | + |r 2 |-2- J R ), (37) 



fi 

which is used together with the WS potential with the " 5 He" parameter set. 
The Hamiltonian employed in our work can thus be written as follows: 

H = (1) + H (2) (38) 

where H^ 1 ' is the one-body Hamiltonian described above augmented by a hard sphere 
Coulomb potential of radius Rq from the 4 He core, and is the two-body interaction 
among valence particles, which can be written as a sum of SGI and Coulomb terms. It 
is important to emphasize that the Coulomb interaction between valence protons can be 
treated as precisely at the shell-model level as the nuclear part. Indeed, the Coulomb two- 
body matrix elements in the HF basis can be calculated using the exterior complex scaling 
as decribed in Ref. 4|. Though not used in this paper, as we are considering only one valence 
proton in the shell model space, this feature of the Gamow Shell Model allows for a precise 
treatment of the Coulomb term. 

The SGI interaction is a compromise between the SDI and the Gaussian interaction. The 
parameter R in Eq. (j3*Tj) is the radius of the WS potential, and V (J, T) is the coupling 
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constant which explicitly depends on the total angular momentum J and the total isospin 
T of the nucleonic pair. A principal advantage of the SGI is that it is finite-range, so no 
energy cutoff is, in principle, needed. Moreover, the surface delta term in (J37j) simplifies the 
calculation of two-body matrix elements, because they can be reduced to one-dimensional 

n 

radial integrals. (In the case of other finite-range interactions, such as the Gogny force [29] , 
the radial integrals are two-dimensional.) Consequently, with SGI, an adjustment of the 
Hamiltonian parameters becomes feasible. 

The Hamiltonian ()38|) is diagonalized in the Berggren basis generated by means of the 
GHF procedure of Sec. IIVI This allows one to use the optimal spherical GHF potential for 
each nucleus studied; hence a more efficient truncation in the space of configurations with a 
different number of particles in the non-resonant continuum. 



B. Choice of the valence space 

The valence space for protons and neutrons consists of the Qp%/2 arid Opi/2 GHF resonant 
states, calculated for each nucleus, and the {^3/2} and {ipi/2} (i = 1, • • • ,n), respectively 
complex and real continua generated by the same potential. These continua extend from 
3?[/c]=0 to 3?[A;]=8fm _1 , and they are discretized with 14 points (i.e., n=14). The Opx/2 state 
is taken into account only if it is bound or very narrow. For the lightest isotopes considered, 
it is a very broad resonant state (r ~ 5 MeV), and, on physics grounds, it is more justified 
to simply take a real {ipi/2} contour, so the completeness relation is still fulfilled. 

Altogether, we have 15 P3/2 and 14 or 15 pi/2 GHF shells in the GSM calculation. The 
imaginary parts of fc-values of the discretized continua are chosen to minimize the error made 
in calculating the imaginary parts of energies of the many-body states. Another continua, 
such as S1/2, d 5 / 2 , ■ ■ ■ are neglected, as they can be chosen to be real and would only induce 
a renormalization of the two-body interaction. We have checked DO that their influence 
on the binding energy of light helium isotopes is negligible. On the other hand, the lsi/2 
anti-bound neutron s.p. state is important in the heaviest Li isotopes ( 10 Li, 11 Li) and plays 
a significant role in explaining the halo ground-state (g.s.) configuration of n Li P,!^. At 
present, however, solving a GSM problem for n Li in the full psd GHF space is not possible 
within a reasonable computing time. 

Having defined a discretized GHF basis, we construct the many-body Slater determinants 



18 



TABLE II: Squared amplitudes of different configurations in the ground state of Li without trun- 
cation (second column) and with truncation to at most two particles in the non-resonant continuum 
(third column). The sum of squared amplitudes of all Slater determinants with n particles (protons 
[it] or/and neutrons [u]) in the non-resonant continuum is denoted by L+. 



Configuration 


No Truncation 


Truncation 


PP3/2M Op§/ 2 M 


0.561-i2.783-10" 4 
0.096+i4.732-l(T 5 
0.184+il.203-l(T 4 
0.064+i2.419-l(T 5 
0.088+i7.032-l(T 5 
0.008+il.62M0 _B 


0.612-i2.285-10" 4 
0.096+i4.980-l(T 5 
0.164+il.077-10" 4 
0.054+il.600-l(T 5 
0.075+i5.508-l(T 5 




from all s.p. basis states (resonant and scattering), keeping only those with at most two 
particles in the non-resonant continuum. Indeed, according to our tests, as the two-body 
Hamiltonian is diagonalized in its optimal GHF basis, the weight of configurations involving 
more than two particles in the continuum is usually quite small, and they are neglected in 
the following. To illustrate this point, let us consider the J n =3/2~ ground state of 7 Li. 
Since, in this case, there are only three valence particle, the complete calculation is possible. 
As seen in Table |Hj the weight of the component with three particles in the non-resonant 
continuum, L^[ttu 2 ] is indeed two orders of magnitude smaller than other configurations 
weights. The presence of pn/ 2 ] modifies the weights of other configurations on a very 
ninor way. For the leading configuration, 0^3/2 MOp^M, without any particles in the non- 
resonant continuum, the effect is ~8%. As a consequence, the neglect of the i}+[t[v 2 ] 
component leaves the overall structure of the state unchanged, and one can safely truncate 
the shell model space while slightly renormalizing the two-body residual interaction. 
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To calculate the first 2 + (T = 0) of 6 Li and in order to have the two first 5/2~ of 7 Li, 
nevertheless, it is necessary to take into account the Opi/2 states, even if they are broad. 
Indeed, without these states, the first 2 + (T = 0) of 6 Li does not exist at the level of pole 
approximation and as a consequence is impossible to find. Moreover, there is only one 5/2~ 
for 7 Li at the level of pole approximation without the Opx/2 states. As a consequence, one 
used the resonant 0pi/ 2 states in the basis to calculate these two states, thus deforming the 
P1/2 contours in the complex plane so the resonant Opi/2 states are enclosed. One also took 
twice as more points for these contours in order to further reduce discretization effects. To 
check that the two set of basis states are equivalent, one calculated the first 1 + , + and 
2 + (T = 1) with them. The first 3 + is not mentioned as no pi/ 2 state, resonant or scattering, 
can enter its decomposition for obvious geometrical reasons. The comparison between the 
two is shown in Fie. where it is clear that the two wavs to calculate the eieenstates are 
equivalent. 

6 . 
Li 



without Opl/2 withOpl/2 



FIG. 5: Eigenstates of 6 Li calculated with and without the 0p\/2 resonant state in the basis. 
C. The lithium chain 

As a pilot example of GSM calculations in the space of proton and neutron states, we have 
chosen to investigate the Li chain. The continuum effects are very important in these nuclei, 
both in their ground states and in excited states. The nucleus 11 Li is also a well-known 
example of a two- neutron halo. In our p-space(s) calculation, we consider the one-body 
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Coulomb potential of the 4 He core, which is given by a uniformly charged sphere having 
the radius of the WS potential. It turns out that the inclusion of the one-body Coulomb 
potential modifies the GHF basis in lithium isotopes as compared to the helium isotopes, 
an effect which is usually neglected in the standard SM calculations. 

Recent studies of the binding energy systematics in the s<i-shell nuclei using the Shell 
Model Embedded in the Continuum (SMEC) (31I Q have reported a significant reduction 
of the neutron-proton T=0 interaction with res pec t to the neutron-neutron T—l interaction 
in the nuclei close to the neutron drip line P. I3 3}] . In SMEC, this reduction is associated 
with a decrease in the one-neutron emission threshold when approaching the neutron drip 
line, i.e., it is a genuine continuum coupling effect. The detailed studies in fluorine isotopes 
have shown that the reduction of the T=0 neutron-proton interaction cannot be corrected 
by any adjustment of the monopole components of the effective Hamiltonian. To account for 
this effect in the standard SM, one would need to introduce a particle-number-dependence 
of the T=0 monopole terms. Interestingly, it has recently been suggested that a linear 
reduction of T=0 two-body monopole terms is expected if one incorporates three-body 
interactions into the two-body framework of a standard SM. 

Our GSM studies of lithium isotopes indicate that the reduction of T=0 neutron-proton 
interaction with increasing neutron number is essential. For example, if one uses the 
Vq(J, T = 0) strength adjusted to 6 Li to calculate 7 Li, the g.s. of 7 Li becomes overbound 
by 13 MeV, and the situation becomes even worse for heavier Li isotopes. To reduce this 
disastrous tendency, in the first approximation, we have used a linear dependence of T=0 
couplings on the number of valence neutrons N n : 

V (J = 1, T = 0) = am, [1 - P w (N n - 1)] , (39) 
V (J = 3, T = 0) = a 30 [1 - Pso(N n - 1)] , (40) 

with aw = -600 MeV fm 3 , /3 10 = -50 MeV fm 3 , a 30 = -625 MeV fm 3 , and 630 = -100 
MeV fm 3 . This linear dependence is probably oversimplified, as shown in Refs. |5|,|33| where 
the proton-neutron T=0 interaction first decreases fast with increasing neutron number and 
then saturates for weakly bound systems near the neutron drip line. For the T=l interaction, 
we have taken parameters Vo(J — 0, T — 1) and Vq(J = 2, T = 1) determined for the He 



ground states 



as they provide reasonable results in the He chain. 



The results of our GSM calculations for the neutron-rich Li isotopes are shown in Fig. |13 
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FIG. 6: Experimental (EXP) and predicted (GSM) binding energies and spectra of lithium isotopes 
obtained with the SGI Hamiltonian. The resonance widths are indicated by shading. The energies 



are given with respect to the core of He. Experimental data are taken from 



aalaa 



One obtains a reasonable description of the g.s. energies of lithium isotopes relative to the 
g.s. energy of 4 He, but excited states are reproduced roughly. Clearly, the particle-number 
dependence of the matrix elements has to be further investigated in order to achieve the 
detailed description of the data. The absence of an anti-bound s\/2 state in the Berggren 
basis is also likely responsible for large deviations with the data seen for 10 Li and n Li. 

In a number of cases, excited GSM states are calculated to lie above several decay thresh- 
olds, i.e., they are predicted to be unstable to, single-nucleon, deutron, proton+neutron 
emission, and/or a decay. The total decay width of a nucleus in a given GSM eigenstate 
is given by an imaginary part of its complex eigenenergy. Different open decay channels 
contribute incoherently to the total decay width and their respective partial width cannot 
be separated easily 2J. In practical applications, however, one may calculate the spectro- 
scopic factors for the separation of nucleon(s) or nucleon groups in a given GSM eigenstate, 
following the well-known procedures of the standard shell model jjjsl, B|, i.e., by calculating 
the probability to find a certain one- or many-particle configuration formed by A — k and 
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TABLE III: Average isospin (|41[) calculated in GSM for various states in 6 ' Li. 



r 


1+ 2+ 2+ 3/2- (g.s.) 1/2- 5/2- 7/2" 3/2~ 


T 

- 1 av 


0.01 0.023+i0.020 1+10.015 0.509 0.514 0.507 0.505 1.503+i0.011 



nucleons in a state of the A-nucleon system. We intend to implement this option in the 
future. 

Due to the explicit presence of the Coulomb potential in the GSM, isospin is no longer 
conserved. In order to assess the isospin-mixing effect, for each shell model state \^) we 
define the average isospin quantum number T av in the following way: 



-1 + a/1 + 4(^|T 2 |^) 
T„ = V 2 (41) 

where the isospin raising and lowering operators in T 2 act on single-proton and single-neutron 
states with explicitly different asymptotic behavior. As seen in Table the values of T av 
indicate very small isospin-mixing effects. The imaginary part for 2 + of 6 Li and 3/2~ of 7 Li 
comes from the fact these states are unbound. This result demonstrates that despite the 
presence of the Coulomb potential and despite a large coupling to the continuum in some 
cases, the isospin quantum number is still almost nearly conserved. That is, isospin is a 
very good characteristic of nuclear states, even if they are unbound, such as the second 3/2~ 
state of 7 Li (which is a T=3/2 isobaric analogue of the 7 He ground state) or the lowest T—l, 
jtt_2+ s tate of 6 Li (which can be viewed as a T=l isobaric analogue of the first excited 
state of 6 He). 

The identification of states can be problematic when the eigenstates have same angular 
momenta and parities. When their isospin is different, the calculation of the approximate 
isospin quantum number is useful to identify one state with its experimental counterpart, 
as was done with the two 2 + states of 6 Li. For the two 5/2" of 7 Li, which have both 
T = 1/2 experimentally, the consideration of their width can be used instead to identify 
them. Indeed, one of them is broad with a width of 880 keV, while the other is narrower 
with a width of 89 keV. As one 5/2" state has to occupy 0pi/2 broad resonant states at the 
level of pole approximation while the other does not, the latter state can be associated with 
the 5/2~ with a width of 89 keV, while the former can be associated with the 5/2~ with a 
width of 880 keV. 
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VI. CONCLUSIONS 

The Gamow Shell Model, which has been introduced only very recently , has proven 
to be a reliable tool for the microscopic description of weakly bound and unbound nuclear 
states. In the He isotopes, GSM, with either SDI or SGI interactions, was able to describe 
fairly well binding energy patterns and low-energy spectroscopy, in particular the Borromean 
features in the chain 4_8 He 0,0]. Using the finite-range SGI interaction made it possible to 
perform calculations in the GHF basis, thus designing the optimal Berggren basis for each 
nucleus. 

In the Li isotopes, the results crucially depend on the T=0 interaction channel. It was 
found that the T=0 force should contain a pronounced density (particle-number) dependence 
which originates from the coupling to the continuum and leads to an effective renormalization 
of the neutron-proton coupling. This effect cannot be absorbed by the modification of 
T=0 monopole terms in the standard SM framework. The effective renormalization of 
(J = 1,T = 0) and (J — 3, T — 0) couplings, and, to a lesser extent, other coupling 
constants found in the present GSM studies has to be further investigated. To better take 
into account all these effects, calculations with a finite-range, density-dependent interaction 
inspired by the Gogny force are now in progress. The three main problems related to 
the realistic GSM calculations are the treatment of the center of mass (essential, especially 



in the context of halo nuclei), the inclusion of anti-bound states |4l|, and the 

handling of very large shell-model spaces. 

For that matter, the successful application of GSM to heavier nuclei is ultimately related 
to the progress in optimization of the GSM basis, related to the inclusion of the non-resonant 
continuum configurations. A promi sing development is the adaptation of the density ma- 
trix renormalization group method (45 ] to the genuinely non-hermitian SM problem in the 
complex-fc plane using the j-scheme 6, l^ . 

In summary, in this paper we report the proof-of-principle proton-neutron calculations 
using the Gamow Shell Model. Our single-particle proton Hamiltonian contains the Coulomb 
term that explicitly breaks isospin symmetry. In order to extend GSM calculations to open- 
proton systems, the Berggren completeness relation has been extended to the case of the 
Coulomb potential. The completeness relation has also been derived for nonlocal interactions 
that naturally appear in the GHF method, a Hartree-Fock inspired procedure to optimize the 
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s.p. basis. According to GSM, the isospin mixing-effects are very weak, even for high-lying 
unbound states of Li isotopes. 
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